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CN ' Abstract 

. In our recent work, we reported an exhaustive study on the simulated bit error rate 

\ (BER) performance of a low-complexity likelihood ascent search (LAS) algorithm for 

detection in large multiple-input multiple-output (MIMO) systems with large number 



ON 



of antennas that achieve high spectral efficiencies. Though the algorithm was shown 
to achieve increasingly closer to near maximum-likelihood (ML) performance through 
simulations, no BER analysis was reported. Here, we extend our work on LAS and 
report an asymptotic BER analysis of the LAS algorithm in the large system limit, 
where Nt,N r — > oo with Nt = N r , where Nt and N r are the number of transmit 
and receive antennas. We prove that the error performance of the LAS detector in 
V-BLAST with 4-QAM in i.i.d. Rayleigh fading converges to that of the ML detector 



^ . as Nt, N r — » oo. 

. Keywords — High spectral efficiencies, large-MIMO detection, likelihood ascent search. 
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Multiple-input multiple-output (MIMO) systems that employ large number of transmit and 
receive antennas can offer very high spectral efficiencies of the order of tens to hundreds of 
bps/Hz [I], [2]. Achieving near-optimal signal detection at low complexities in such large- 
dimension systems has been a challenge. In our recent works, we have shown that certain 
algorithms from machine learning/ artificial intelligence achieve near-optimal performance in 
large-MIMO systems that employ tens of transmit and receive antennas using V-BLAST 
and non-orthogonal space-time block codes (STBC) [3] with tens to hundreds of dimensions 
in space and time, at low complexities. Such algorithms include local neighborhood search 
based algorithms like a likelihood ascent search (LAS) algorithm [I], [5] and a reactive tabu 
search (RTS) algorithm [5J, and algorithms based on probabilistic data association (PDA) [7] 
and belief propagation (BP) [E],[S]- Similar algorithms have been earlier reported in the con- 
text of multiuser detection [TUJ-JTB]. In through detailed simulations, we have shown 
that LAS and RTS algorithms achieve increasingly closer to maximum-likelihood (ML) per- 
formance and that PDA and BP algorithms achieve near maximum a posteriori probability 
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(MAP) performance for increasing number of dimensions in large-MIMO systems. For e.g., 
in [5], the BER performance of the basic LAS algorithm (which uses a single symbol update 
based neighborhood definition) and its generalized version (which uses a multiple symbol 
update based neighborhood definition) has been exhaustively studied through simulations. 
However, BER performance analysis of the LAS algorithm for large-MIMO detection has 
not been reported. In this correspondence, we fill some of this gap by presenting an asymp- 
totic BER analysis of the LAS algorithm in the large system limit, where N t , N r — > oo with 
N t = N r , where N t and N r denote the number of transmit and receive antennas, respectively. 
Asymptotic performance analysis of large systems in the context of multiuser detection and 
MIMO communication have been reported in the literature [I7]-[22|, using random matrix 
theory (e.g., PZ],[IE|), replica method (e.g., [Tj5] , [2D] , [2T] ) , and free probability theory (e.g., 
[22]). We, in this correspondence, present an asymptotic BER analysis of the LAS algorithm 
in the large system limit. Specifically, we present an analytical proof that the error perfor- 
mance of the LAS detector for V-BLAST with 4-QAM in i.i.d. Rayleigh fading converges 
to that of the ML detector as Nt, N r — > oo with Nt = N r , which is an analytical result that 
has not been reported so far. 

The rest of this correspondence is organized as follows. The MIMO system model and 
LAS detection algorithm are summarized in Section [2j The asymptotic analysis of the LAS 
algorithm is presented in Section [3j Lengthy proofs of lemmas and theorems are moved to 
the appendices. Simulation results and discussions are presented in Section |U Conclusions 
are given in Section 

2 System Model 

Consider a V-BLAST erystem with Nt transmit antennas and N r receive antennas, Nt < N r . 
Let x c £ C NtXl denotqjthe symbol vector transmitted, and H c £ C NrXNt denote the channel 
matrix such that its (i,j)th entry hij is the complex channel gain from the jth transmit 
antenna to the ith receive antenna. Assuming rich scattering, we model the entries of H c 
as i.i.d. CM(0, 1). Let y c £ QN r xi anc j n ^ e ^N r xi denote the received signal vector and 
the noise vector, respectively, at the receiver, where the entries of n c are modeled as i.i.d 



1 Vectors are denoted by boldface lowercase letters, and matrices are denoted by boldface uppercase letters. 
[.] T and [.] H denote transpose and conjugate transpose operations, respectively. ||.|| denotes Euclidean 
distance. 
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CJ\f(0, a 2 ). The received signal vector can then be written as 

y c = H c x c + n c . (1) 
Let y c , H c , x c , and n c be decomposed into real and imaginary parts as follows: 

y c = yi + JYQ, x c = x 7 +jxQ, n c = n I +jn Q , H c = H 7 +jH Q . (2) 
Further, we define H r . e R 2 ^ x2Af ', y r e R 2NrXl , x r e R 2NtXl , and n r £ M 2 ^ xl as 

H r = ( h q "h? ) ' Yr = [yJ y 2 T) Xr = [x ' X 5 ]T> nr = [nJ n 2 T (3) 

Now, (Q]) can be written as 

y r = H r x r + n r . (4) 

Henceforth, we shall work with the real- valued signal model of the system in (jlj). For 
notational simplicity, we drop subscripts r in (jlj) and write 

y = Hx + n, (5) 

where H = H r G M. 2N - x2Nt , y = y r e E 2JV - xl , x = x r G M 27VtXl , n = n r G M 2Ar - xl . In this 
real-valued system model, the real-part of the complex data symbols will be mapped to 
[x±,--- ,XN t ] and the imaginary-part of these symbols will be mapped to [xN t +i,--- ,x2N t ]- 
For M-QAM, [xi,--- ,XN t ] can De viewed to be from an underlying M-PAM signal set and 
so is [a;jv 4 +i , ■ ■ ■ ) X2Nt\- Let Aj denote the M-PAM signal set from which Xi takes values, 
i = 1, 2, • • • , 2N t ; e.g., for 4- QAM, A; = {1, -1} for i = 1, 2, • • • , 2N t . Now, define a 2N t - 
dimensional signal space S to be the Cartesian product of Ai to A 2 jv t . The ML solution 
vector, cIml, is given by 

argmin 2 arg min , T T T . 

= de§ ll y ~ Hd " = deS ( dHHd - 2 y Hd )- ( 6 ) 

In the following subsection, we summarize the low-complexity LAS algorithm, using a neigh- 
borhood definition based on 1-symbol updates, presented in [5] for large-MIMO detection 
for M-QAM. The channel matrix H is assumed to be known perfectly at the receiver. 

2.1 LAS Algorithm for Large-MIMO Detection 

The LAS algorithm starts with an initial vector d^°\ given by d^ ' = By, where B is the 
initial solution filter, which can be a matched filter (MF) or zero-forcing (ZF) filter or MMSE 



filter. The index m in d^ m ' denotes the iteration number in a given search stage. The ML 
cost function after the kth iteration in a given search stage is given by 

C« = dW T H T HdW-2y T Hd( fe ). (7) 

The d vector is updated from fcth to (k + l)th iteration by updating one symbol, say, the 
pth symbol, as 

d (fc+D = d « + AWe pi (8) 

where e p denotes the unit vector with its pth entry only as one, and all other entries as 
zero. Since d( fc ) and d^ fc+1 ^ should belong to S, Ap can take only certain integer val- 
ues. For example, for 16-QAM, A p = {—3,-1,1,3}), and X p can take values only from 
{—6, —4, —2, 0, 2, 4, 6}. Using ([7]) and (jHJ), and defining a matrix G as 

G = H T H, (9) 

we can write the cost difference C( fc+1 ) — as 

jr(/<*>) t cW-CW = l^a p -2l^\z p % 

where z^f* is the pth entry of the vector given by = H T (y — Hd( fc )), a p = (G) pp 
is the (p,p)ih entry of the G matrix, and l p k ^ = \X P \. The value of l p which gives the 
largest descent in the cost function from the /cth to the (k + l)th iteration (when symbol p 
is updated) is obtained as 

„(*) 



/(*) 



2 a p 



(10) 



where [•] denotes the rounding operation. If dp were updated using l p k lp t , it is possible that 
the updated value does not belong to A p . To avoid this, we adjust l p k l pt so that the updated 



(k) 

value of d p belongs to A p . Let 



arg mm ^,iu\ , , , 

p 



If o t ) < 0, the update for the (k + l)th iteration is 

d (*+D = d« + ^sgn(^))e s (12) 
z (*+i) = z ( fc )-/i5 pt sgn(^))g s , (13) 

^s,opt) 



where g s is the sth column of G. If ^{l^lpt) — 0, then the search terminates, and is 



declared as the detected data vector. 
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3 Asymptotic Analysis of LAS Algorithm 

In this section, we prove the asymptotic convergence of the error probability of the LAS 
detector to that of the ML detector for N t , N r — > oo with N t = N r in V-BLAST. Consider 
4- QAM, i.e., § £ {+1,-1} 2A \ and let N t = N r . An n-symbol update on a data vector 
d £ § transforms d to (d — Ad n ) such that (d — Ad n ) £ S. Further, (d — Ad n ) is obtained 
by changing n symbols in d at distinct indices given by the n-tuple u n = ^2, ■ • • ,%), 
1 < ij < 2Nt,Vj = 1, • • • , n and ij ^ ik for j ^ k. Therefore, we can write Ad n as 

n 

Ad n = ^2d ifc e ijk , (14) 

k=l 

where d{ h is the i^th element of d. Let L n C § denote the set of data vectors such that for 
any d £ L n , if a n-symbol update is performed on d resulting in a vector (d — Ad n ), then 
||y — H(d — Adn) 1 1 > ||y — Hd||. Our main result in this section is Theorem [2j To prove 
Theorem [21 we need the following Lemmas [I] to [51 Slutsky's theorem [23], and Theorem [TJ 
Lemma 1 Let d £ S. Then, d £ L n if and only if, for any n-update on d, n £ [1, 2 ■ • ■ , 2N t ] , 

(y-Hd+^HAd n ) T (HAd n ) > 0. (15) 

Proof: By definition, if d £ L n , then no n-symbol update can result in a reduction in the 
ML cost function. Using this, we can write 

||y-H(d-Ad n )|| 2 > ||y-Hd|| 2 . (16) 

Simplifying ( 1161) . we get ( Tl5i) . Since the choice of the indices in u n is arbitrary, the lemma 
holds true for all possible n-tuples of distinct indices. For the converse, if d satisfies f TT5"j) for 
all possible u n for a given n, then, since ( fl5l) and ( fl6l) are equivalent, d also satisfies ( fl6l) for 
all possible u n . This implies that d £ L n . □ 
If d £ Li , then using Lemma [T] and §5§ , we can write 

(n + H(x - d) + h p d p ) T (h p d p ) > 0, Vp = 1, ■ • ■ , 2N t , (17) 
where h p is the pth column of H. 

Lemma 2 Assuming uniqueness of the ML vector &ml in (0|) ; a symbol vector d £ § is the 
ML vector if and only if the noise vector n satisfies the following set of equations 

n j, n 

(n + H(x-d)+(j>A)) (5>A)>0> (18) 
i=i 3=1 
Vn = 1, • • • , 2N t , and for all possible n-tuples (ii, ■ ■ ■ ,i n ) for each n. 
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Proof: If d is the unique ML vector, then from the definition of the ML criterion in ([6]), it 
must be true that any n-update on d will not result in any decrease in the ML cost function. 
Therefore, d G L n , V n — 1, 2 • ■ ■ , 2N t . Hence, by Lemma [1], it must be true that d satisfies 
f fT5|) for all n — 1, 2, ■ ■ ■ , 2N t and for all possible u n for each n. Substituting y = Hx + n in 
f fl5|) . we get (fT8l) . This proves the direct result. To prove the converse, let the noise vector n 
satisfy ffl8|) for some vector d. Since y = Hx + n, the conditions in ffl8l) imply the conditions 
in ( fl5l) for all n — 1, 2, • ■ • , 2N t and for all possible u n for each n. Therefore, by Lemma (TJ 
d G L n for all n — 1, 2, • • • , 2N t , which then implies that d indeed is the ML vector. □ 

Definition: For each d G S and for each integer m, 1 < m < 2N t , we associate the set of 
vectors TZ dm = {v | v G R m and (v + H(x - d) + ( £J =1 h v^)) T (E?=i h *A) > 0, Vn = 
1, • • • , m, and for all possible n-tuples (ii, • ■ ■ , i n ) for eachnj, and define 7?.d = ^d 2iV t • 

Lemma 3 7/ i/ie noise vector n G 7?.d; then d zs i/ie ML vector. Let dj, dj G § and dj 7^ dj. 
Then TZ di and TZd^ are disjoint. 

Proof: From Lemma [2] and the definition of TZd, it is clear that d is the ML vector if and 
only if n G TZd- The disjointness of TZ^ and Jt^, i 7^ j, can be shown by contradiction. If 
TZdi and are not disjoint, then there exists some vector v belonging to both 7Zd t and 
TZdj • If v were to be the noise vector n, then, v would satisfy the set of equations in (fTBj) for 
both d = dj and d = dj, since v belongs to both 7Zd t and TZd 3 , This, by Lemma [2], implies 
that both dj and dj are ML vectors, which is a contradiction because of the uniqueness of 
the ML vector. □ 

Lemma 4 Let h G M. 2Nt be a random vector with i.i.d entries distributed as A/"(0,0.5). 
Let {hj},z = 1, 2, • • • , m be a set of vectors, with each hj G M? Nt and having i.i.d entries 
distributed as W(0, 0.5), E[hjhJ] = for i 7^ j , and E[hhJ] = for j — 1, • • • , m. Then 

l im ^k=i n n fc = (19) 

N t ->co mN t 

Proof: Let h = ^= Y^k=i hfc- Then, h ~ jV(0, |). Therefore, we have 

lim E ^ hThfc = lim J£*L. (20) 

iv t ^oo mN t N t ^oo ymN t 



We can write 



lim ^ = lim S£^, (21) 



where hk and hf. are the kth elements of h and h, respectively. The r.v's h^h^ k = 
1, ■ • • , 2N t are i.i.d with mean zero. From the strong law of large numbers [23], it follows 
that lim^^ Y? k Z\ z% = °- Using this in ([20]) completes the proof. □ 

Before we present the next lemma, we present the Slutsky's theorem on convergence of 
random variables, which is used to prove Lemma [5] and Theorem [TJ 

Slutsky's Theorem J2$j : Let {X m } and {Y m } be sequences of random variables. If {X m } 
converges in distribution to a random variable X, and {Y m } converges in probability to a 
constant c, then it is true that i) {X m + Y m } converges in distribution to X + c, ii) {X m Y m } 
converges in distribution to cX, and Hi) converges in distribution to 

Lemma 5 For a given u n and a given d £ §, define a r.v z Unj d cls 



where ij £ u n , j — 1, • • ■ , n. For any u n and any d G §, z Un ,d converges to zero in probability 
as N t — > oo, i.e., z Un ,d — — ^ as N t — > cxd, = 2, 3, • ■ ■ , 2iV t . 

Proof: Proof of this Lemma is given in Appendix A. □ 

In Fig. [T] we plot the simulated pdf of -2 Un ,d for n = 2N t for different values of N t = N r for 
a certain u n and d (the pdf was observed to be same for different u n and d). We observe 
that with increasing N t = N r , the pdf of -2 u „,d tends towards the Dirac delta function at 
zero. This implies that z Un ^ tends to zero in distribution, and hence in probability, for large 
Nt — N r , which is formally proved in Lemma [51 

Theorem 1 Let d £ S and n £ TZ^i. Then n £ TZ^ in probability as N t — > oo, i.e., for any 
8, < 8 < 1, there exists an integer N(5) such that for N t > N(8), p(n £ TZ^) > 1 — 5. 

Proof: Proof of this theorem is given in Appendix B. □ 

Theorem 2 The data vector/bit error probability of the LAS detector converges to that of 
the ML detector as N t , N r — > oo with N t = N r . 

Proof: Let diAS be the final output symbol vector of the LAS algorithm given x, H and n. 
The algorithm terminates if and only if no 1-update results in any further decrease of the 



2, 



u n ,d 



E 




(22) 
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cost function. This implies that for the given x, H and n, d^ A s £ L 1; and therefore it must 
be true that n satisfies ( TTT|) with d replaced by d^^s- These set of equations are the same 
which define the region TZ^i. Therefore, replacing d by diAs, we can equivalently claim 
that n G TZ d i^ A . Using Theorem [H we can further claim that asymptotically as N t — > oo, 
n G 7^d i4S in probability From Lemma [31 we know that if n G lZd LAS , then d^s is indeed 
the ML vector for the given x, H and n. Therefore, we can state that asymptotically as 
Nt — > oo, diAs is indeed the ML vector in probability. That is, for any 5, < 8 < 1, there 
exists an integer N(5) such that for N t > N(5) 

P(d LAS is the ML vector) > (1 - 5). (23) 

Therefore, we can write that for N t > N(5) 

PLAs(error) = P(d LAS ^ x) = P{d LAS ^ x | d LAS = ML vector)P(d LA s = ML vector) 

+ P(d LAS ± x | d LAS ^ ML vector)P(d LA5 ^ ML vector). (24) 

From (1231) . we have P(di A s ^ ML vector) < 5. Also, P(dL A s ^ x | d^^ = ML vector) is 
the probability of error for the ML detector, which we denote by PML(error). Using these, 
we can bound the probability of error for the LAS detector as 

PLAs{zrror) < P ML (error) + 5 P(d LAS ^ x | d LAS ^ ML vector) < P ML (error) + 5. (25) 

Since 5 can be arbitrarily small, we can conclude from (125]) that indeed as Nt — > oo, the 
symbol vector error probability of the LAS detector converges to that of the ML detector. 
This proof can be adapted to show that apart from the symbol vector error probability, the 
bit error probability of the LAS detector also converges to that of the ML detector. The 
proof for the bit error probability convergence is along the same lines as ( 1241) and (12"5"|) , except 
that instead of defining the error event as d^ A s ^ x, we define error events for each bit. For 
example, for the pth bit, the error event is defined as d PLAS ^ x p . □ 

4 Simulation Results and Discussions 

In Fig. [2] we show the simulated BER performance of the LAS detector for V-BLAST with 
4-QAM and MMSE initial vector for increasing N t = N r . Since an analytical expression for 
ML performance in the large MIMO system limit is not available and simulating the ML 
performance for large dimensions involves prohibitively high complexity, we plot the SISO 
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AWGN performance as a lower bound for comparison. It can be seen that for increasing N t = 
N r , the BER performance of the LAS detector approaches the SISO AWGN performance 
at high SNRs. Figure [3] shows the average SNR required to achieve a BER of 10 -3 for 
increasing N t = N r and 4-QAM. It can be seen that, for large N t = N r , the required 
SNR gets increasingly closer to that required in SISO AWGN for increasing N t = N r . A 
similar behavior can be observed in Fig. |4] for 16-QAM as well. In Figs. [3] and HI we also 
see that there is an initial degradation in performance for increasing number of antennas 
(for Nt < 10). This shows that the LAS detector is suboptimal for small systems with 
small number of antenna^, and becomes optimal in the large system limit (as proved in the 
previous section). LAS detector achieves close to large system limit performance in systems 
with large number of dimensions (e.g., hundreds of dimensions in Figs. [3] and [3]). Such large 
number of dimensions need not be realized in spatial dimension alone, as in V-BLAST. As 
shown in [5], exploiting time dimension in addition to space dimension, large non-orthogonal 
STBC MIMO systems can render large dimensions with less number of transmit antennas 
that can be implemented in practice. A 16 x 16 non-orthogonal STBC from cyclic division 
algebra [3] with complex data symbols has 512 real dimensions; with 64-QAM and rate-3/4 
turbo code, this STBC achieves a spectral efficiency of 72 bps/Hz. In [5], LAS algorithm has 
been shown to achieve near-capacity performance in 16 x 16 STBC MIMO systems even in 
the presence of spatial correlation and with estimated channel matrix. Further, considering 
that NTT DoCoMo has demonstrated a 12 x 12 V-BLAST MIMO system operating at 5 
Gbps at a spectral efficiency of 50 bps/Hz at 10 Km/hr mobile speeds [23], the availability 
of low-complexity large-MIMO detection algorithms like the LAS algorithm analyzed in this 
correspondence can motivate the adoption of 16 x 16 and 24 x 24 MIMO systems operating at 
spectral efficiencies in excess of 50 bps/Hz in emerging wireless standards like IEEE 802.11 
VHT and IEEE 802.16/LTE-A. 

5 Conclusions 

We conclude with the following two remarks: i) The derivation of analytical BER expres- 
sions for the ML performance in the large MIMO system limit for different signal sets is an 
open problem. Since large MIMO systems can be viable in practice due to the availabil- 

2 We do not have a theoretical explanation for this small system behavior of the LAS detector, whereas 
we are able to prove its asymptotic large system behavior. 
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ity of low-complexity detectors like the LAS detector, analytical BER expressions for the 
ML performance in the large MIMO system limit would be quite useful as a benchmark for 
comparing the performance of practical detectors in large-MIMO systems. The statistical 
mechanics approach employed in [19] for large CDMA system BER analysis can be inves- 
tigated for such an analysis, ii) While we are able to prove the asymptotic convergence of 
LAS performance to ML performance for 4-QAM here, our simulation results for higher order 
QAM (e.g., 16-QAM; see Fig. H|) show similar behavioral trend like that for 4-QAM. Conse- 
quently, we conjecture that such a convergence holds for general M-QAM and an analytical 
proof to show this can be attempted as an extension to this work. 

Appendix A: Proof of Lemma [5 

We present the proof of Lemma [5] in this appendix. The proof is by mathematical induction 
on n. Base Case: For n = 2, we have to show that 



h r h 

dpdq W¥ E +Wp JU ° asJV *- >0 °« V ^9 = l,2,"-,2iV t , p^q. (26) 



h T h 

We can write the random variable ,,, J\ M2 as 

l|h P ||- + ||M- 



h p T V(2iV t ) (27) 



(HM a + IIM 2 )/(2iv*)' 

As Nt — > oo, by strong law of large numbers, the denominator of (1271) converges to 1 almost 
surely. Also, the numerator of (12 7p can be written as 

1L p ll 1 2—ik=l ' l P,ki l q,k 



2N, 2N t 



(2t 



where h p k and h q ^ refer to the /cth entry of the vectors h p and h q , respectively. Each 
h P} khq } k term in the summation in (1281) has the same distribution and has mean 0. Therefore, 

h T h 

by strong law of large numbers, we can see that 9 converges to almost surely. This 

h T h 

also implies that £ N q converges in distribution to the constant 0, and hence by Slutsky's 

h T h 

theorem, ^ [J^pp || 2 conver g es m distribution to 0. Since, if a sequence of r.v's converges in 
distribution to a constant then the sequence converges in probability to that constant, we 

h T h 

conclude that indeed || h ||2 P + ||h |p converges in probability to 0. This proves the the base case. 
Induction Hypothesis: Let z Un ^ — as N t — > oo, Vn = 2, 3, ■ ■ • , m. 
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Induction Step: Proof for n — m + 1: We have 



2 U(m+i),d — y^^ 1 II Vi II 2 ^ ' 

Z^j=l \\ n ij\\ 



|h II 2 + V m lib- II 2 



E^xUh^.H 2 ^ E 3 m =1 l|h ;j IP 

J _|_ H hi (m+l)H 2 



(30) 
(31) 



E7=illN'H 2 

Using Slutsky's theorem and the strong law of large numbers, it can be shown that the 
denominator in ( 13TT) converges to (1 + — ) in probability. Also, from the induction hy- 
pothesis, the term — - — =yn — p^]p — - — in the numerator of ( 13T]) converges in probabil- 
ity to 0. Therefore, the numerator in ( l3Ti) converges to the same distribution that the 

'3=1 l|n V 



term ^„T m converges to. Also, the term ^ H , u2 is the same 

2~,j=l ll ni j II ^3=1 H n, 3 II 

as r^m^nu msT// — itf> • Further, from the strong law of large numbers, the term 

(Sj=i l|h«jl| 2 )/( m -^) converges almost surely to 1. Therefore, from Slutsky's theorem, we 
know that ^m 1 ',,, now/ — kf\ converges in distribution to the distribution to which 

(Ei=i \\ 2 )/{mNt) to 

the term (J2k=i hJ {m+1) h ik d Hm+1) d ik ) / (mN t ) converges. 

For a given vector d, hi k di k is a random vector whose distribution is the same as that of 
h ik . Therefore, applying LemmalU we see that the term (^^ = x^, m+1) ^i k d^ m+ ^di k ) / {mNt) 
converges almost surely to 0. Hence, the numerator in (l3~lj) converges in probability to the 
constant . Therefore, z U(m+1) ,d — —> as N t — »• oo. This proves the induction step and 
completes the proof of Lemma □ 

Appendix B: Proof of Theorem ffl 

We present the proof of Theorem [TJ in this appendix. We shall prove through induction that 
if n £ then n £ TZ^m in probability, Vm = 2, ■ ■ • , 2N t , as N t — > oo. Base Case (m = 2 j: 
Let n £ ^d 1 - Therefore, from the definition of TZ^, n satisfies ( 1T7|) . We show that n £ 7^d 2 
in probability as — > oo. For n to belong to 7?.d 2 ; m addition to satisfying (|17[) . n must 
also satisfy the following equation Vp, q — 1, ■ • ■ 2N t ,p ^ q: 

(n + H(x-d) + h p d p + h q d q ) T (h p d p + h q d q ) > 0, (32) 
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which can be rewritten as 

(n + H(x - d)fh p d p + (n + H(x - d)) T \i q d q > -\\h p \\ 2 - \\h q \\ 2 - 2d p d q h^h q . (33) 

Since n satisfies (jTTl) . it satisfies the following two equations: 

(n + H(x-d)) T hpd p > -||h p || 2 , 

(n + H(x-d)) T h^ > -||hj 2 . (34) 



Comparing (1341) and ( 1331) . we notice that if h p and h q are orthogonal, then n trivially satisfies 
(13"3"|) for all N t . Therefore, when h p and h q are non-orthogonal, the only extra term in the 
RHS of ( 1331) is 2d p d q hph q . Applying Lemma [3, with n = 2, we see that as Nt — * oo, the r.v. 

h T h 

— — converges to zero in probability. Then, we can write, for any e, < e < 1 



IIM 2 +||h 9 

p ( iiJmLii' > f ) K m > /(f) ' (35) 

Now, let us analyze p(n £ TZ d ?) for the case of d p d q = +1 (a similar analysis holds for d p d q = 

— 1). Consider two disjoint events E\ = j [[hjp+jfhgii 2 < e } anc ^ ^ 2 = { j|h TT^Tp p > e }' 
Then, we can write 

p(n£ft d2 ) = p^^^p^) + p(n £ 1l d 2\E 2 ) p(E 2 ) . (36) 

The event i?i can be further split into two disjoint events En and E\ 2 , given by En = 
{0<h^h g <e(||h p || 2 + ||h g || 2 )} and £ 12 = {0 > hjh, > -e (||h p || 2 + ||hj 2 )}. Also, from 
(l35|) . p(Ei) > 1 — e and p(E 2 ) < e. Therefore, using (1361) . we can write 

p(n^TZ d2 ) < p(n^TZ d2 \E l )p(E 1 ) + e 

< p(n £ K d2 \E n )p(En) + p(n £ TZ d2 \E 12 ) p(E 12 ) + e 

< p(n £ ^l^u) + p(n £ ft d2 |£ 12 ) + e. (37) 

If event E n is true, then 

-(||h p || 2 +||hj 2 ) > -(||h p || 2 +||h f/ || 2 + 2hjh f/ ) > -(||h p || 2 + ||hj 2 )(l + 2e). (38) 
Since n £ TZ d i, n satisfies (1341) . and hence satisfies the following equation: 

(n + H(x-d)) T (V P + V ? ) > -(||h p || 2 +||hj 2 ). (39) 
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Using (1381 and fl39|) . we see that n satisfies fl33j) . and therefore n e 7£ d 2. Hence, we can 
conclude that p(n ^ 72.^2 |-En) = 0. Now, we can rewrite fl37|) as 

p(n^^ d2 ) < p(n^ft d2 |£ 12 ) + e. (40) 

If event E12 is true, then 

-(IIM 2 + ||hJ 2 )(l-2e) > -(||h p || 2 +||h ? || 2 + 2hjh 9 ) > -(||h p || 2 + ||hj 2 ). (41) 
Using fl33l) and (jHJ), we can write that 

p(n^ d 2|£ 12 ) = p([(n + H(x-d)) T h p d p + (n + H(x-d)) T h g d g < -||h p || 2 - ||hj 2 - 2hjh 9 ] | E 12 
< p ([ - (||h p || 2 + j|h 9 j| 2 ) < (n + H(x - d)) T (h p d p + Vy < - (||h p || 2 + ||h,|| 2 ) (1 - 2e)\ \ E 12 ). (42) 

Define 1Z € to be a set of vectors in R 2Nt , as 

^ = {v I - (||h p || 2 + j|h f; || 2 ) < (v + H(x - d)) T (h p d p + h q d q ) < - (||h p || 2 + j|h ? j| 2 ) (1 - 2e)} . (43) 
Also, define a function f 2 as 

/ 2 (e) = p(neft e |£ 12 ). (44) 
Using the above definitions, (H2|) can rewritten as 

p(n ^ TZ d 2 I £ 12 ) < / 2 (e). (45) 



Let 61,62 G K, 61,62 > 0, and e\ > 62- From the definition of lZ t in (143]) . it can be seen 
that C 7Z ei . This implies that /2(ei) > f2{ti)- Hence / 2 is a monotonically increasing 
function. Using (1451) . we can rewrite (1401 as written as 

p(n i Tl A 2) < / 2 (e) + e. (46) 

Therefore, 

p(ne^) > l-(/ 2 (e) + e). (47) 

Now define (72(e) = /2(e) + e- So ^ 2 is a monotonic function and is therefore invertible. Let 
5 = (72(e). Using (|35|) and the above definitions, we can write that 

N t > /(e) 

> N 2 (5), (48) 
13 



where N 2 = / o g 2 1 . We can then write ( 1471) as 

p(neK d 2) > 1-8. (49) 

Since g 2 is a continuous monotonic function, for any 8, < 8 < 1, there exists an integer 
N 2 (5) such that for N t > N 2 (8), p(n £ 7?. d 2) > 1 — 5. Therefore, n £ 7Z d 2 in probability as 
A" 4 — > 00, thus proving the base case. 

Induction Hypothesis: Let n £ T^™- 1 in probability as N t — >oo. 

Induction Step: We need to prove that n £ TZ d m in probability as Nt — > 00. For n to belong 
to IZ dm , n must satisfy the following equation for all possible m-tuples (ii, i 2 , ■ ■ • , i m ): 

(n + H(x-d)+(5]h i .d i .)) (X>Aj > 0, (50) 
i=i ^ i=i ' 

which can be written as 

m—l 

(n + H(x - d)) r ( M%) + ( n + H ( x " d))% m d im 
i=i 

> -IIESx'H^II 2 - l|h im || 2 -2(Er= 1 1 h i ,d i ,)Vd im . (51) 

However, we know from the induction hypothesis that (n + H(x — d)) T (^™7 1 1 h^d^) > 
~\\ Y^j=i hijdijW 2 . Also, since n £ U d i, we know that (n + H(x - d)) T h im d im > -||h im || 2 . 
Therefore, if the term 2(^™i~ 1 \\i j di j ) T h.i rn di m in the RHS of (ISTp were 0, then (1501 would have 
been trivially satisfied. We now show that the contribution of the term 2(^™~ l 1 h^di^h^d^ 
when compared to the other two terms in the RHS (15TT) converges to as N t — > 00. 



Define a r.v. t> m = S-i v, ' J . 'ii nT T2 • Our objective is to show that as A^ — > 00, v. 

II Ej=l "-ijdijW + ll"i m l 



24.Hl,. 112 " WUJ - ^"-"J^^"^ ^ ouajvv vllujv ' i/ m 



A || E"^l h i rf i IP 

in probability. This is equivalent to proving that w m = v m + 1 — 13 %J — 



INJI 2 +UE™T S^H ; 
converges to one in probability as N t — > 00. We can write W m clS 



"P m llh- II 2 
E-=i llhijIP + E^ill^H 2 

From LemmaEl we know that for any integer m, 1 < m < 2N t , it is true that — M . 3 m — - — 

E/=l ll n ijll 

converges to in probability as A^ — > 00. By Slutsky's theorem, this is equivalent to 



2 spm y^-m uTu J J II spj=m i > || 2 

+ 1 = v4n n r no — ► 1 53 



V m llh- II 2 V m llh 

Z^j=l ll^jll ^.7 = 1 H 
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as N t — > oo. We shall use this result to prove the convergence of w m in (152]) . Using (1531) . it 
can be seen that the numerator of w m in (1521) converges to 1 as N t — >• oo, i.e., 

" Sr^^f 1, asiV^oo. (54) 

In the denominator of (1521) . it can be shown that the term 

llh,J| 2 p 1 



as iVj — > oo. (55) 



/j=i n i 

The 2nd term in the denominator of (1521) can be rewritten as 



z^=m .,n E-r/iih^ir 



(56) 



Similar to the derivation of ( |53l) . we can claim that the numerator in (|56l) converges to one 

llh' II 2 1 

in probability. From Slutsky's theorem, it can be shown that |^=?jjh~"jj5 conver S es ^° i R 
probability. Using this and Slutsky's theorem, it can be shown that fl56l) converges to 
in probability. Using this result along with (p>4"j) . fl55l) and Slutsky's theorem in fl52l) . it can 
be shown that w m converges to one in probability as N t — > oo. This, therefore, implies that 
v m converges to zero in probability. As proved in the base case, it can be shown that for any 
5, < S < 1, there exists an integer N m (5) such that for N t > N m (5), p(n e TZd^) > 1 — 5. 
This proves the induction step and completes the proof of Theorem [TJ □ 
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Figure 1: Simulated pdf of z Ull ,d for n = 2N t for increasing N t = N r . 4-QAM. The pdf tends 
towards Dirac delta function at zero. 




- A 


- Nt = Nr = 1 


- e 


-Nt = Nr=10 




-Nt= Nr = 50 


- + 


-Nt= Nr= 100 


- $ 


-Nt= Nr = 200 


- B 


-Nt = Nr = 600 




- SISO AWGN 



7.5 



8.5 9 9.5 10 

Average Received SNR (dB) 



10.5 



Figure 2: Simulated BER performance of the LAS detector for V-BLAST as a function of 
average received SNR for increasing values of Nt = N r . MMSE initial vector, 4-QAM. LAS 
detector achieves near SISO AWGN performance at high SNRs for large Nt = N r . 
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Figure 3: Average received SNR required to achieve a target BER of 10~ 3 in V-BLAST for 
increasing values of Nt = N r for 4-QAM. LAS detector with MMSE initial vector. LAS 
detector achieves near SISO AWGN performance for large Nt = N r . 
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Figure 4: Average received SNR required to achieve a target BER of 10 -4 in V-BLAST for 
increasing values of N t = N r for 16-QAM. LAS detector with MMSE initial vector. LAS 
detector performance approaches SISO AWGN performance for large N t = N r . 
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